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SPIRAL MOTIONS OF VISCOUS FLUIDS* 


By Georg Hamel 


INTRODUCTION 


The equations for the plane motion of viscous fluids of constant 
volume are, after elimination of the pressure and introduction of the 
| stream function V by which the velocity components 
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are expressed, reduced to the one equation 


SAV. DAV AV OAV _ ! 
eo ae oe ee oe (T) 


therein o indicates the ratio between viscosity coefficient and spe- 
cific mass u, and A signifies the Laplace operator. 


This equation is satisfied by all potential motions 


| Ay = 0 


Ih (4 

: however, this fact is of little significance since viscous fluids adhere 
i to solid walls and, from well-known considerations of function theory, 

| there cannot exist a potential motion which would do so. Otherwise, 

Ä properly speaking, only Poiseuille's laminar motion is known as exact 
solution of equation (I) and that solution does not even show the sig- 

3 nificance of the quadratic terms because they identically disappear 
there. 


1 *"coiralformige Bewegungen zäher Flüssigkeiten."  Jahresber. d. 
| deutschen Math. Ver. 25, 1917, pp. 34-60. 
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Under these circumstances it seems perhaps useful to know a few 
more exact solutions of equation (I) for which the quadratic terms do 
not disappear; such solutions will be indicated below according to two 
methods. 


In both cases, one deals with motions in spiral-shaped streamlines 
(which are observed frequently). 


Third, we shall, in addition, investigate the neighborhood solu- 
tions to pure radial flow. 


FIRST PART 


We raise the question: 


Are there solutions of equation (I) which are not potential motions 
for which, however, the stream paths are the same as for a potential 
motion whereas the velocity distribution is to be different? 


We shall be able to indicate such solutions, in fact all of them: 
the streamlines are logarithmic spirals (including concentric circles 
and pure radial flow); for the velocity distribution, one arrives at an 
ordinary differential equation which for pure radial flow leads to 
elliptic functions. In the discussion, the influence of the quadratic 
terms becomes manifest in a considerable difference between inflow and 
outflow (see paragraphs T, 8, and 9). 


We require, therefore, solutions y of equation (I) for which 


y = f(ọ) 


and AP = O, but not Ay = 0. The latter condition excludes 


f''(Ọ)=0 
We limit ourselves to steady motions 2r = Oy 


1. The calculation becomes clearer if first the auxiliary problem 
has been solved: 
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Transformation of equation (I) into isometric coordinates, that is, 
t such curvilinear coordinates %9,%X that 


| p + ix = w(x + iy) = w(z) 
| Let us thus assume 

| 

N 

| 

| 


= oP _ ox dp. ox 
v(o, x) Sx > 


2 2 
If one denotes SV + SV by A'Y, there results first, with the abbre- 


| ap? 0X 
= viation 
à 
dw| -= 
dz a 
Av = Q A! Y 


_ With the double integral extended over an arbitrary region, one has 


J) (3 33 - ZA hax ay - J| aaa 
- [Ren ar ALAN lap ax 


XGA’) oF. aa) 27 òP ox _ AP OX 
- J = 3% ot Sy % > ax d 


Since, however, 
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there follows 


SAW _ IAS _ g (ne oy 34A av) + a (3.12.8 ov _ Ding a 
Ox Oy Oy ox cp ox ox OP op ox ox op. 
However, 
dw 
9 In Q. 5 Og in We oR Lin W- on 4 
op op d dw dw V? 
(az) 
and 


2 
dw 
Jing -2.9 Rin Æ- -2 È Jin- -27 Din u -2g -4 
Ox X dz o d dw d (=) 
Z 


are valid. If one puts the analytic function of z 


one obtains 


QAQU QOQAOYV. 


dw 
dz 


n 
SA! SV OA OV if, ov ov 
(A Sx OX at) + ara Bev DD) 
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Finally, there results 


CUN AU 


op X ox 


N 


Bid: oe; y A'S o AY - BA'Y v) 
dla AV + AV E + (9:2 a - LA »)| 


e 
ln Q = ln En is & harmonic function, thus 
A' In Q@=0 
hence, 


arg = (Bina), (Dana) = a? + 22 
Q op OX 


Thus one obtains as the result of the conversion of equation (I) to 
isometric coordinates %9,%X for steady motion 


APP 
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Per 


A aav 09 - 2a 34 + arya +» BY) =o A AV + 
i SD. 3X ax Ə ax 9 


Am " 
E XOT ae, 
= 


la? « v2 SAY SAY 
i A'y(a +2) + oe Sx j 


| therein, a + bi is the analytic function 


and A' denotes the operator 


a a 


QAN(QA'Y) = A A'T + QA'QA'Y + aaa N 9 A'V = 


Li 


(11) 
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2. We return to the question On page 2: y must be a mere function 
of € 


y = £(P) 


If derivatives, with respect to Ọ, are denoted by primes, equation (II) 
becomes 


f''f'b- o| 27V + £'' (a + b) + apina | (IIT) 


f may depend only on ® but must not depend on X. 


This is certainly possible if a and b do not depend on X, thus, 
since a + bi is an analytic function of ® + xi, do not depend on 9 
either, if a + bi is, therefore 


that is, constant. We shall see later (paragraph 3) that this is the 
only possibility. 


From a and b being constant, there follows 


W =- —— In (z - 20) + Wo 


thus, after introduction of the polar coordinates 


Z2-2. = revi 


p santa ln r + b8) +9 
2 2 O 
a +b 


NACA T 1342 2 


Thus, the streamlines = const are identical with the logarithmic 
spirals 


SS en 


a lin r + bð = const 


a = 0 signifies pure radial flow, b = O flow in concentric circles. 
The velocity đistribution, however, is given by equation (III): the 
radial component is 


_ av r: OP _ 28a f: L 
or Or 5 eee r 
consequently pude LP l the magnitude of the velocity. There- 
2 2 
a +b 


fore, f' must disappear on solid walls. 


Without restriction of the generality, one may presuppose left- 
hand spirals so that r increases with %, thus a and b have dif- 
ferent signs; since, furthermore, -( + iX) is an analytical function 
just as Ọ + iX, and equation (III) is actually invariant with respect 


1 


to a simultaneous sign” change of P,a,b, one may presuppose 


oJ 
HA 
O 


Therefore, positive velocity components ov and 8 signify for 


DU outflow,: 
in contrast for 


f' «0 inflow. 


t 


LTranslator's note: The original says "time change,’ 
misprint. 


obviously & 


8 NACA TM 1342 


Since € may be replaced by cf, one may in addition impose a con- 
dition on the constants a and b. 


3. We now want to conduct the proof that on the basis of our require- 
ments a and b must be constant, that therefore the flows in loga- 
rithmic spirals are the only ones the flow patterns of which correspond 
to a potential motion without themselves being a potential motion. 


If a and b were not constant, the analytical function a + bi 
would produce a conformal transformation of the Ọ + iX-plane; by virtue 
of equation (III) which with the abbreviations 


(S Bu; LV 
hui B=4 f' Cafi 
» E 20 Tu 
(f'' 20 is excluded) may also be written 
a“ + b - 2A(P)a - 2B(P)b + C(D) = O (III') 


the circles (equation (III')) would correspond to the straight lines 
Ọ = const in this transformation. 


These circles would therefore have to form an isometric curve 
family. 


However, if the family of curves 


(III') 


I 
© 


g(a,b,?) 


is to be an isometric one so that AY = O, the function g must satisfy 
the equation 


^gg, _ 2&, (Ep aba + Zy xp) + & o (Ba = e?) B” (IV) 


fo = 25+ 2 
dat db“ 


and this equation must either be identically satisfied, or be a conse- 
quence of equation (III'). 
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One has 
| 


~ 6, =2(a-A), & - 2(b - B), Ag 


h, Bpa = -2A', Ry = -2B' 


a 


Bp = -2A'a - 2B'b + C', Ep = -2A''a - 2B''b + C'' 
g^ * g^ = (a - A)? + h(b - B)? = h(A? + B? - C) 


Thus, equation (IV) is quadratic in a and b; an easy calculation shows 
the result that the quadratic terms are automatically eliminated. There- 
fore, the coefficients of the two terms must be zero whence follow three 
conditions 


_A''t 0! - 2AA' - 2BB' _ B'' _ Ct - 2AA' - 2BB' 
1 H 
A "P a B DN RR 
_c'' | œ! - 2AA' - 2BB' 
t 
C g- A -p 


n Hence, there follows further that A', B', C', must be proportional and 
; furthermore that 
| 


ii B! 
N Qus AC aoe 
i 
i: must be constant. 
n 
if The final condition yields 
i 
i C -aB* A=y7,B - 5 
| or with 
a y 
“L zz Q fd = =y 
20 20 


: pIV - ati! + pri! and Fer = YI LS + of! 
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which, integrated, ylelds 


TIE: = l af'? + BL! + € and prr = ER yf'? 4 of! + n 


D9 
tà 


Comparison of the two values for f''' results in 
= of? + pf£' +e 
fre = 2 
yf' + 8 


which must be identical with the preceding value of f''. The compari- 
son requires 


y =0 a =0 B = 5° € = 5n 


thus C = A? = 8^ constant and f'' = Bf' + n. 


The second condition 


NEN ER = Dr. = const 
à à à 


A +B -C B 


however would yield — = const and this together with f'' = Bf! +n 
x b 
would result in the contradiction 


f' = const 


therewith, the proof has been produced. 


4. We now turn to the determination of the velocity according to 
differential equation (III) which may be integrated once and assumes, 
after introduction of the quantity proportional to the velocity at unit 
distance 


u = f'(p) 


P Ax - — 
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—— 
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the form 
u'' + 2au! + u(a? + v2) - D uê +C =0 

20 


This equation is identical with a damped oscillation which takes place 
under the influence of the potential 


-B u3 +4 2 (a? + bê) u + Cu 
60 2 
We start with the limiting cases: 
l. The streamlines are concentric circles: b= 0. 


Then 


and, because of 


u = const + r(A + B4 ln r) 


whereby, the velocity distribution 


is given. The exact solution of Conette's case is also contained therein: 
the three constants here are determined from the two limiting values of 
the velocity and from the fact that in case of a full turn around the 
circular annulus, the pressure must revert to its initial value. An 

easy calculation yields B4 = O and thus 


v = S(r? -~ r3?) 
r 1 


(More details on the determination of the pressure are seen in para- 
graph 10.) 
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2. Ihe flow is purely radial: a= Q. 


The differential equation reads 


and leads to elliptic functions 


u: = up t 30bu? + const u + const 
E l- 2. (e1 - u) (ez - u)(e, - u) 


where the three e's are only subject to the one condition 


ei + e t e. = 30b 


but otherwise are still at disposal. 


Since, according to the remark on page 8 one relation between a,b 
is still unused, it will be expedient to put 


b= -2 
so that one obtains, according to page T 
9-238 
Then the conditional equation for the e reads 
ej + e% + ej = -60 


and one has 
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thus 


= u = -20 +P) ——( - 30); E2783 


Vés 


where 89585585 are the three integration constants. For the pressure 


(see paragraph 10) there results the equation 


SÍB,.1,9N.1 gp + 20 gii I Lfl pi2 , 208) 
opu 2 > = | 09 „2\2 


its uniqueness is & priori ensured, thus does not determine here any of 
the constants. 
Discussion of the Radial Flow 


o. The condition 


e, + e, + e, = -60 (1) 


requires at least one e to have a negatively real part, for instance 


R(e,) 2 R(e,) 2 R(e3) 


then 


with the equality sign being valid only when all three e's have the 
same real part. 


Furthermore, since this part is real, there must apply 


(a) for three real e's either 


-o < u Í e £ -20 
3 
or 


eo £ u £ e 
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(b) for one real e 


where, however, this e may be positive. 
Furthermore, two possible types of flow must be distinguished: 


l. Either there are no solid walls, thus a source or sink in an 
unlimited fluid. Then u must be a periodic function of 9, with a 
period which is an integral part of Zr. u = -œ is excluded, u= 0 
need not occur. Therefore, this case can occur only for three real e's, 
and 


must be valid. 


2. Or there are two solid walls, for instance for 3 = O and for 
81 (which may also be equal 2n); then at these walls u must be 


il 


u = 0. 


(a) In case of three real e's there must be, additionally 


< > 
en £ O e g O 
and either 
(a) PEA 
or 
(B) O<u< e4 


(b) In the case of one real e, this e must be positive and 


One remembers, furthermore, that according to page 7, paragraph 2, 
u>0 signifies outflow, u <0 signifies inflow; so that one has 
inflow in the case of 2(a)(a), and outflow in the case of 2(a)(ß), and 
2(b) above. For the case 1, both cases may occur. 
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First Case: Free Flow 


6. One must assume 9=0 for u= e, and has therefore 


9 =. B9 du 
in a RIDE 


Hence, there must 


[I mmm 2x (2) 
e jo n 


2 ler c 9 - 93 7 93) 


with n being an integral. 


By the known substitution 


e -e 
u=e, + (e - e ) sin“ X^ = Ł £ 
2 l 2 en -e 
2 3 
equation (2) becomes 
x 
2 is av . [£z 
3o n 


E O 
vo too Vi + X*sin^v 


If one now introduces the mean velocity? 


2At the distance r= l. 
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there becomes because of equation (1) 


m -4 
ep - 8, = 60 + zu, SED (3) 
thus 
X^ = 5 - 
60 + Zum - 5 ô 
and 


x | l 
ET dV nn 60 + au, - 55 (2') 
n 

z l + X?sin*V V 9c 


From this, one may draw several interesting conclusions. 


One has 
14 x x 
I ay |. du ay "ET ay 
e O O 
Vi + X*sin*w Va 4 x sin y l + X*cos^ 
uj 
" Ih ——— 
O 1.2 
1+ 5% (1 - cos 2Y) 
x 
l >a | dy 
1.2 0 1.2 
1+ 5% (1 + cos ZW) 1+ 5H 
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thus 
atone - x 

2 [2 ay 2 x 

x 

O 
l + X*sin^VW Ni +S ex 
where 
Ocgcecl 


Thus the relation (2') between up 95 Dn,0 reads, due to the significance 


of x^, 
— Óá—— a 
l n 60 
\ [6a + 3u4 - zu - €)ö 
or 


2 
60 + Zum - = n8 = 60 T (2'') 


with n being a proper fraction. 


Since, furthermore 


x x 
I x ay - I? x av = sinn E 
O O 
1 + X^sin^V i eae 


thus becomes arbitrarily large with increasing X, one has lime = O, 
2( 200 


thus, 
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From equation (2'') there follows 


2 
n 
Un > aep - = 


which, with u = l, gives as the minimum value 


m>- o 


The mean inflow velocity is therefore considerably limited upward, the 
more so, the easier movable the fluid. 


However, this is the only restriction: If U, and the integer n 


are selected so that 


n 


2 
En 


there exists, certainly, a pertaining Ö. 


For if "E increases from zero to the value 60 + 3u,, = n*^85 


lies between zero and 60 + Zum (because for the second value x 


becomes infinite and, hence, n° = 1) so that certainly sometime = n^b 


becomes equal to 60 + 3u, - 60 T which is presupposed to be positive. 


One sees, furthermore, that for a prescribed period number n the 
fluctuation 5 and for a prescribed fluctuation 6 the period number n 
must increase to infinity with the mean velocity u,- 


Second Case: Outflow Between Solid Walls 


7. The cases 2(a)(ß), and 2(b) may be summarized thus 


jg = E a Sn 
2 O 2 
VC u u) (u + 2au + B) 
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e >O is the maximum velocity (at the distance r = 1); because of 


2a = -e, - e, = 60 +e 


and ß = ee, >O, otherwise, however, arbitrary 


u^ + Zau + B 


may for prescribed e assume all values from ue 


e 
vy = 2 | ee eee | RET 
V e 
O 


(e - u) (v4 + Zau + B) 


+ Zau to o, so that 


appears not at all restricted downward, but upward restricted by 


I = ? "e — 
l,max V2 
O \/(e - u)u(u + e + 60) 


since 


one has 


9 = 2n SEN PENA (4) 
l,max 2e(l + e) + 120 


where e signifies a positive proper fraction. 


For the outflow, the width of the wall opening appears therefore 
restricted, according to the preceding equation, by the maximum value e 
of the velocity. For small velocity and large viscosity, the maximum 
lies near n, otherwise, however, lower; with increasing e it drops 
below all limits. 
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If, therefore, an angle opening smaller than x is prescribed, it 
permits an outflow only up to a certain maximum value. If a greater 
outflow quantity is prescribed, the jet will, therefore, actually prob- 
ably separate from the walls. 


Also, there is, of course, for any prescribed angle a flow 


9 
1 
possible where partly inflow partly outflow occurs. 


Third Case: Inflow Between Solid Walls 


8. There remains the case 2(a)(a) 


e, Suso 
all three roots e real, e3 e negative, e positive. 
Here 
O 
z 39 du 
94 = 2 7 
e " = = 
2 (a en) (u €3)(*1 u) 
O 
= [és | du 
E 
(u- e -u^ - 2au - p) 
2 
where 


otherwise, however, arbitrary. Thus, the angle 94 may be made arbi- 


On the other hand, however, it may 


the 


trarily small for prescribed en. 
also be made arbitrarily large: one takes, for prescribed ZE 
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negative value e sufficiently close to en, 88 far as this is not made 


3 
impossible by ei < 0. The sole relation between the e 


e. te. t e. = -6g 


however, results with ei 2-0 in 


-e. - e_ > 60 


à 3 


If -€ Z 30, e, may actually be assumed arbitrarily close to en. 


3 


If the maximum inflow velocity is larger than 30, any angle 94 


is possible between the solid walls. 
If, however, “En < 30, say €30, where € is a positive proper 


fraction, only 


-e4 =e, + (6 - 3e)o = -e, te, + 6(1 - e)o 


is possible and 


= Y G 
(u e?)|» * (6 3e)o + e; |(eı u) 


attains its highest value for ej = O 


O 


31 ‚max = Vee zu 


i NC + 3eco)|u + (6 - 3e)o| C2) 
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where n and 5 are positive proper fractions. Thus, the maximum of 
94 is larger than m. 


When the maximum inflow velocity is smaller than 30, the angle 
openings of the solid walls also may attain any magnitude up to s. 
Flow in Spir&ls 


9. Because of the damping 2au' (see paragraph 4, page 11), a 
periodic solution, aside from u = const, is not possible. 


A free motion in logarithmic spirals is always a potential motion. 
In contrast, there exist other flows on logarithmic spirals between 
solid walls. 

In order to have, for r = const, the variable (p agree with the 


angle ð, one may furthermore prescribe for the constants a,b, in such 
a manner that one obtains 


P + b^ 


thus 


a must be a proper fraction, otherwise it remains arbitrary. 


Equation III, once integrated, (see page 10) then reads 


2 
u'' + 2au' + 8^u + E u? +C =0 


2 F an} sa <4, but Sa. 


z 
+ 
(b 
H 
(D 
To 
l 
| 
D 
o 
l 


rr 
"t: 
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The velocity at unit distance is 


If one first omits the damping, one has exactly the same case one 


had before except that 
2 
B instead of 2: 
60 3c 


is in front of the square root (see page 12). The relation for the e 


remains the former one. Since g^ < 4, the angle opening is increased 
by this influence 94 - 


The damping, however, takes effect in the same sense. Nevertheless, 
the main result remains correct. 


For outflow the admissible angle opening 9 is restricted by the 


23 
maximum flow velocity in such & manner that it tends toward zero when 
this velocity increases beyond &ll limits. 


If one puts 


-a(p 


u= ve 


the above differential equation becomes 
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If © =O is assumed to be the location of the maximum Yo for v, 


multiplication by 2v' and integration yields 


2 v 


m 
vts (p? = a” )(v" - v 2| +È vee Pay + 2C e av = 0 
z ao Jy v 
l O 


From the corresponding equation for u 


u 2 
12 ^ (u^ d P (u? 3) 
+ ha + = + = + = - 
u I u' du + f^ (u uo gzu ug ec( uo) O 
O 
one can see that for equal uo the u-curve will be the steeper, 94 
therefore the smaller, the larger C. For value close to ug this is 


immediately clear from the differential equation, for in case of u' = 0, 
u'' will be the smaller, the larger C, thus |u'| the larger; from 
the preceding equation one may see, however, that for larger C 


u 


u uo 
u'? 4 ha | u'duc-u'?- s f u' du 
"0 


has the higher value. Hence, follows directly for the ascending branch 
(u' > 0) that always |u'| has the higher value when C is the larger 


value. For if u' would once reach for the initially flatter curve 


u 
(C smaller) the value of the steeper curve, [ u' du would have to 
uo 
40 
have for the former the smaller, thus L u' du the higher value 
u 


which for u' > O immediately leads to contradiction since, up to then, 
that is between u and uo; u' had been the smaller value. If, how- 


ever, u' < O, one has in case of a variation of the C by A 


u 
A —U_-- + he i | ? Au! du = 2 AC 
u 


— 


SS A Se eee 


_— al 2 
eU UELUT. 


NACA TM 1342 25 


or, since 


du'® = AJu'|(21u'| + Alu’ |) 


' uo 
alu (ale) t Apul), a f A|u'| an + 2 AC 
bo De u. - U u 


If there were at one point Alu'| = O, then at this point the first 
term would, for fixed AC, decrease with decreasing u, that is, go over 
from positive to negative values; the second term also would decrease 

u 
since the part of the integral f > Alu'| du supervening with 
u 


decreasing u would be negative. This is impossible, however, since 
the sum of both terms is supposed to be constant 2 X. 


Therewith, it has been generally proved that the angle opening 94 


decreases with increasing C (for fixed uo); since the maximum possible 


9 is desired, the minimum admissible value for C may be assumed. 


l 


This value is determined from vie > 0 
MA 2 v 
2C 9 geras 2 - (e - a?) (v A v?) -BL S vve ag 
v » 20 v 


whence, one may see that the minimum admissible value of C is zero or 


negative? ‘ 


The above inequality must be valid for all v's between Vo and 
zero and for the positive and negative P attained. One may write them 


2 -a(P, + P 
ddo - (8° " a”) (vo x weer " Wo” + vov + v^ e "(C2 o) 
g 


` 


P xscordinE to page 22, pe - a? >0. 
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where PoP are certain mean values. One must further note that, for 
Y = Vos [P| and 2 | mist be zero whereas they have maximum values 


for v - 0. The severest restriction is due to the absolutely smaller 
value of the right side, thus the minimum admissible C is given by 


= ! = t t 
2C - - (8° - a^)v e “To - p v fe 5% ucl 
0 6g 9 
Po" Po. are positive and the maximum of Po and P, which occur for 
y = 0. 


Consequently, one has for the maximum possible 4 


-aP,' pO 
via = (e? - a?) (vo? - v?) - Voe 0 ] e&av + 
v 


2 YO -2(9,'+9,' vo 
B [ vee aay - ACE : o') H eS ay 
v 
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and since 9,' < 84 


whence follows 


(Compare formula (h), page 19). 


Hence follows that with increasing Yop thus also with increasing 
Ups ð must drop below all limits: a certain width of the spiral 
permits only a limited outflow velocity. 


SECOND PART 


lO. Sinee the only spiral motion, possible without walls, of the 
type used so far, lead to be a potential motion, exact steady and non- 
steady two-dimensional motions in free spirals will be investigated 
according to another method. 


In polar coordinates the differential equation (I) reads 


9A, 1/9 Al OY _ ð Ay O9 = oA AV 
ot or af oP or 


where 


> 
N 
hi jr 


Drd pl 
ror 9 rê y 
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Obviously this equation permits solutions which are linear in € 
vV=ut Ox 


in order to make the velocity which has the components 


vy 2X and V -U ç 9X 
r j or Ar 


unique and thus enable a free motion, # must be constant. By this 
statement, the differential equation becomes 


Oo Au 
ot 


+ 


H [X 

A od 
E 
! 
Q 
D 
E 


with 


Here it is necessary also to investigate the pressure lest perhaps 
in a motion about the singular point r =0 a multivaluedness of the 
pressure becomes evident. 


Now one may write the equations of motion without elimination of 
the pressure in the following form 


Or because of the invariance of the last term 


3 (9 A p (s A as EM 
= AV dU + A t dr - Aea g ap 


r oP or 
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Hence, follows 


| S[B ,1.2 - ay t. 2( - è) -= _ m 3 Au eu 
(= + v2) NS r o A St X Au - or = un 


By virtue of the differential equation for u the right side is con- 
stant; thus it must be zero to make the pressure in case of a revolution 


i about r= DO revert to its former value, so that one obtains 
à 
r ð Au = .9u. .X Au = O 
or 9 or ot 
which by introduction of r au - y assumes the form 
r 


oy. (a + ah oe = Be (V) 


Steady Motions 


11. The solution independent of t is 


4 when X = -2, otherwise, when X = b 
uy o o 


E 2 
u = c, (In r) +c, lnr +c, 
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If one disregards the trivial case of potential motion, a spiral motion 
the velocity of which disappears at infinity exists when 


A+1<o 
G 


that is, X < - o, thus a sufficiently strong inflow takes place. 


The spirals then have the form 


If &+2< O, they approach at infinity the logarithmic spirals; near 
g 


the sink, in contrast, they converge considerably less pronouncedly 
toward the sink point, and the vortex velocity is considerably higher 
than in case of potential flow in logarithmic spirals. 


Unsteady Motions 


12. If one uses the formulation 
nt 
v-e X, (r) 


one Obtains from equation (V) 


ct 
- 
wi 
Ui 
Se 
£ 
fe 
ct 
- 
ct 
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o 
o 
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Oc 
R 
(D 
< 
D 
ct 
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where the J are the Bessel functions 


a 
J ~2y)\= const r|14.2I 2. 0d 4 
^ c 4 1(1 +2) 


SE) —— — +. . 
gf 2'(1 +A)(2 + X) 


If X does not happen to be an integer, rs, and ry, may be 


regarded as independent solutions. 


13. Similarly to the case of the heat conduction equation there 
exist also of equation (V) integrals which show for r=0 and t = 0 
an indeterminate point. 


Since the differential equation (V) remains unchanged if v is 
multiplied by an arbitrary factor, r by a similar factor, t by its 
square, there must exist solutions of the form 


2 
v=-r £) = rH By (2 
=u (= At Pw (z) 


After substitution, one obtains for w the differential equation 


wit + wat A) + uat - 2a 4 B| L0 (VI) 
- Mz - 


When does this equation permit & solution of the form 


w = eM? 


A simple calculation yields 
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and then either 


or 


a = 2X B-2--X-1 
One thus has two simple integrals of the required type 


Ir 
hg t 
y = tà le 


and 


l rê 


r 
ho t 
y = pl. 


for X = O, both are transformed into the known integral of the heat 
conduction equation. 
Let us continue the discussion of the differential equation (VI). 


The singular point z = 0 is a determinate point. The determining 
equation reads 


2 
p? + pla- X) + em = 0 


and has the roots 


© 
i 
> 
bo Jg 
o 
Il 
i 
NIR 


nen 
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so that generally there exist developments of the form 


Xr e 
2 2 
W Z =Z, l+c_z2z+ce2”"” +... 
1 ( T 2 ) 
and 
„& 
E 2 1 T 
Wo = % (2 + Cy zZ + ens zf +... | 
thus 
œ 
B-A+ = 2 4 
v = r^ 211 + Cy —~— +e r RA ur 
hot 2 
(lot) 
and 
a 
B+= 2 4 
y. t “1+ c, ' m + c, r 7 t. : 
(Lot ) 


with the power series continuously converging since z= O is the only 
singular point of the differential equation. 


If one assumes fp = O, that is, if one desires solution of equa- 
tion (V) of the form 


rê 
iu Gs) 


an integration by definite integrals is possible. 


The differential equation (VI) reads after introduction of the 


roots 
P4 ^9. 


2 
ada Ww dw " _ = t 
er + zZ 2 Q PL” Po + z) + WP} Po O (VI') 
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The connection with Gauss’ equation for the hypergeometric function can 
be easily recognized. If one makes the Euler transformation 


W= feh - Sy y(s) de 


with the integral extended over a suitable closed path, one finds for y 
a differential equation which may be satisfied for 


-l+p, 
n = -P) by y=s8 
and for 
-ltpi 
n= -p by y=s 
à 
Therefore 
Pi -1+0 
w= et «d S ds 
Z 
and 


are integrals of equation (VI'). The integrals are extended best over 
a path which leads from R(s) = +% around the points s - 0 and s=z 
back to R(s) = +o. 


Since 


-p, -l*p 
fei - s) ts Cas 


RR — NEN | 


u ee 


= - = 


ae qe eee ee 


Uk 


DNE, GC ee 


7 
— un 


en — — 77 
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+ ae 0D 
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is analytically regular in the neighborhood of z = O, there is 


-P1 -1-«p 
zi -3/1 - E a 
Wy je ( E S ds 


and 
-l+p 
- -3/1 . & l 
Wo c, [e ( B S ds 


One can show th&t 


e» eel a (c) ) «ce s) 


a= = œ 


are the general solutions of equation (V) and likewise are represented 
by definite integrals in closed form. I shall perhaps refer back to 
this and to the connection with the representation and the development 
in terms of Bessel functions elsewhere. 
THIRD PART 
Neighborhood Solutions to Radial Flow 


14. We shall first look for steady neighborhood solutions to the 
radial flow (pages 12 and 13) by putting 


= f(0) + p(P,r) 


where p is assumed to be a small quantity, the square of which is 
neglected. 


We then obtain for f the former equation 


" IV) 4 agi: er: 
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with f' =u and integrating once 


u''+ho+tu2+¢=0 
o 


For p we obtain 


where 


Since the differential equation in the steady case remains unchanged 
when r and p each are multiplied by an arbitrary factor, there must 
exist solutions of the form 


p = r^w(9) 


One obtains for w the differential equation 


wlV + Ve -W+h+2u- à u) + Zu w' + 
g o 


2 3 
vib = af Bam yy ba) = 0 (VII) 


We are particularly interested in free flows and thus in periodic solu- 
tions in €. 


As concerns the uniqueness of the pressure (see paragraph 10, 
page 28), one obtains by & simple calculation, the condition 


o pax 
“al wla - (a - 2)o]ap = o 
O 
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On the other hand, there follows from the &bove differential equation 
itself, by integration over the interval from O to 2x with assump- 
tion of periodicity 


er - 2) |" wlu -(x- 2)o| ap = 0 
O 


so that in general the uniqueness of the pressure follows from the peri- 
odicity except for the case when à = 2. 


For free flow, u itself is a periodic function of €; the period 
is an integral w part of 2x. However, it is not necessary that w 
have the same period as u; but this period must likewise be an integral 
part of 2x. 


Since 


as well as 


u'? MS - hue - > u3+D= Ze - u)(u - e) (u - e4) 


may be rationally expressed by u, it will be useful to introduce u 
instead of ® as independent variable in equation (VII). Because of 


wt 29W qs 
du 
2 
wow. + GW qe: 
à à 
witt W u' + 3 IH yy + Wyre 
du du en 
CENT, 3 a a“ d 
wiv = d wW oq: +6 IH yt4yrt 4 39 u'' 4+} W u'u''' + QW „IV 
4 3 2 2 du 
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and because of 
uV = (^ - Eun E ee (A a2 Ju? 
o 0 oO 


all coefficients of the new equation are integral and rational in u; 
indicating the degree, one writes them 


4 3 à 

Ro SH + Re SH + Ry IH HR, Ù + Row = 0 (VII') 
à à 3 2 du 
u du du 


with 


Re - u'# = 6 - u)” (u - e3) (u - ej)" 


From the form (VII) one can see that w possesses singularities 
only where they occur for u, thus certainly not in the real part of 9 
(which is of interest); equation (VII') shows that, as a function of u, 
w becomes singular only at the branch points ti e 


Since Rs = 6u'“u'' is divisible by (ei - u)(u - e,)(u - G2); 
the points 1292593 are determinate points, and since the degree of 


the coefficients decreases steadily by 1 with the order of the deriva- 
tives, u = æ glso is a determinate point; the differential equa- 
tion (VII') belongs to the Fuchs class. 


A well-known calculation yields as the four roots of the determining 
equation for the points e the values 


- 3 
2 


Although, therefore, two root differences here are integral, no loga- 
rithmics appear ir the developments: For from the form (VII) there fol- 
lows that at the points Y for which u becomes = e, where, therefore, 
u and u' are regular functions of 9, w also must be such a func- 
Lion, whereas ln(u - e) does not possess this regular character. 
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Therefore, the solutions of equation (VII') have at every point e the 
form 


u = P (u - e) + VE - eP,(u - e) 
other singularities do not exist in a finite domain. l 


For u = œ there results the determining equation 
2 242 3 
2 = 
( u^ + u 3) (4 u x) O 


- which has the roots 


I E H =} 


which are independent of X, and the roots 


which are dependent on i. 


Continuation 


| 

15. Solutions with the real period 2x (this period must be present 
at least in case of free flow) will exist only for certain A. In 
analogy with Hermite's method for Laine's differential equation, one can 
proceed as follows: 


If Wy» 


w(® + 2x) are expressed homogeneously linearly by the w 


WWW, are a fundamental system of equation (VII), the 


h 
v, (9 * 2x) gn ay WP) (v = 1,2,3,4) 
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There certainly exist periodic functions of the second kind, that 
is, there exist solutions w for which 


w(® + 2x) = aw(P) 


This a is a root of the equation of the fourth degree 


a i 7 8e 535 "d E 
a21 "aeg T On Bes E 
D(a; A) = = 0 
“31 832 ck - (x ey 
"M he 243 "d o9 


If a periodic solution is to exist, a = 1 must be a root, and one 
obtains for X the equation 


D(1,3) = 0 


The characteristic exponents which were calculated suggest the 
attempts of putting 


w=u+ const, W= Jn -e and wW = (e, 7 Y) - es) 


Elementary calculation yields the following particular solutions: 
l. The trivial possibility w-u for X-0 


2. w-u for X = 2, that is 


y = f(9) + or*f' (o) 
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where a must be small &nd therefore with the same approximation 


y = £ (9 + ar“) 


80 that the streamlınes are approximately the spırals 


f' remains the same elliptic function discussed before in the case of 
radial flow. It as true that this flow now cannot exist as free flow, 
since this is precisely the exceptional case à = 2 (see page 37); and 
the condition for the uniqueness of the pressure can certainly not be 
satisfied for w= u. 


3. w - u * 30, when X -1 and C = 30 whence for e, - e, < o3 


no contradiction results. 


h. w= ju- e, for X= -1 and e, = O. This solution has a 


period twice that of u; likewise, w = Je -u for à = -l and e = 0. 


Oo. we (2, 7 yu - e3) or We Ju - €5)(B - e3) when A =1 


and e, = O or e = O. This solution too has a period twice that of u. 


The large X may be easily calculated approximately from equa- 
tion (VII). For such large X there is in first approximation 


a à 


IV w'' -Aw-0 


w t 2X 
that is, w= et? (we restrict ourselves to the periodic solutions), 
go that ex is the period. The large set-apart -values are therefore 


approximately integral. 


Finally, one case may be calculated quite elementarily: the case 
when u is constant, the basic flow therefore an all around uniformly 
distributed flow. 


42 NACA TM 1342 


This case is also of significance for the more general one since 


according to a well-known theorem by Cauchy and Bobtznuna the period 
of w in first approximation is obtained if the constant mean value is 
inserted for the periodic u, under the presupposition that the larger 


fluctuation ei - €, be sufficiently small. 


For constant u there follows from equation (VII), page 36 


WIV 4 wit (an? = a +h £u) twit = a3 t aË t a2 £u) = 


thus with the formulation 


la ar) ext a3 t aË ta Bd umo 
g 8] 


This equation has four roots 
u= dr w= (X- 2) ~-a2Ay 
(ej 


so that all integral positive and negative AX are possible (potential 
motions) as well as all A which are calculated from 


with integral uu. 
For the case u = const, that is: 


16. For the radial flow which is uniform all around, the unsteady 
neighborhood solutions also can be given. 


"Boltzmann, Ges. Abh., Bd. 1, S. 43. 
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43 
The differential equation now reads (see page 36) 


or Sere 


One may integrate it either by means of the formulation 


Ap = ext mit r) 


(n integral) and thus arrives at the differential equation 


which may be solved by Bessel functions, or by means of the formulation 
(compare page 32) 


; 2 . 
- enl® r _ ani 
Ap =e (= eli rp ( z) 


whereby one obtains for w(z) the differential equation 


h " hoc] s 


2 2 
ze! ulm + TID ee at > 
o 


For 


zZz =0 this equation has a determinate point, the determining equa- 
tion has the real roots 
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By introduction of the roots and Po the differential equation 


Pi 
assumes the form 


This is, however, exactly the differential equation (VI') of page 33 so 
that everything said about it there is also valid here. 
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